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This chapter is devoted to a discussion of quantum phase transitions in 
regularly alternating spin-^ Ising chain in a transverse field. After re- 
calling some generally-known topics of the classical (temperature-driven) 
phase transition theory and some basic concepts of the quantum phase 
transition theory I pass to the statistical mechanics calculations for a 
one-dimensional spin-^ Ising model in a transverse field, which is the 
simplest possible system exhibiting the continuous quantum phase tran- 
sition. The essential tool for these calculations is the Jordan- Wigner 
fermionization. The latter technique being completed by the contin- 
ued fraction approach permits to obtain analytically the thermodynamic 
quantities for a "slightly complicated" model in which the intersite ex- 
change interactions and on-site fields vary regularly along a chain. Rig- 
orous analytical results for the ground-state and thermodynamic quan- 
tities, as well as exact numerical data for the spin correlations computed 
for long chains (up to a few thousand sites) demonstrate how the reg- 
ularly alternating bonds/fields effect the quantum phase transition. I 
discuss in detail the case of period 2, swiftly sketch the case of period 3 
and finally summarize emphasizing the effects of periodically modulated 
Hamiltonian parameters on quantum phase transitions in the transverse 
Ising chain and in some related models. 



1. Classical and Quantum Phase Transitions 

In this chapter I consider the effects of regular alternation of the Hamil- 
tonian parameters on the quantum phase transition inherent in a one- 
dimensional spin-i Ising model in a transverse field. Before starting to dis- 
cuss this issue let me recall some common wisdoms from statistical physics. 
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2 O. Derzhko 

One of the aims of statistical physics is to describe different phases which 
may occur in many-particle systems and the transitions between different 
phases, i.e. the phase transitions. We often face phase transitions in every- 
day life. Melting of ice, boiling of water or vanishing of the magnetic proper- 
ties of iron after heating are well-known phenomena for everyone. Normally 
we associate the changes in properties of a substance as the temperature 
varies. However, such changes may occur at a fixed finite temperature while 
some other parameter (such as presssure) varies. 

P. Ehrenfest proposed a classification of phase transitions. In the Ehren- 
fest classification of phase transitions we say that the phase transition is of 
the order one if the free energy is continuous across the phase transition 
whereas its first derivatives with respect to temperature and other vari- 
ables (for example, pressure) are discontinuous. Similarly, we say that the 
phase transition is of the order two (three) if the free energy and its first 
derivatives (first two derivatives) are continuous across the phase transition 
whereas the second (third) derivatives are discontinuous. A few years later 
L. Landau proposed to consider discontinuous or continuous phase transi- 
tions. In the Landau sense the discontinuous (continuous) phase transition 
is characterized by a discontinuous (continuous) change in the order param- 
eter and, therefore, it is usually viewed as of the first- (second- or higher-) 
order. In the first-order phase transitions the two phases coexist at the phase 
transition temperature. Thus, if ice is melting one observes two phases, i.e., 
liquid and solid, which coexist at the temperature of phase transition. This 
is an example of discontinuous phase transition. In the second-order (and 
higher-order) phase transitions the two phases do not coexist. An exam- 
ple is the Curie point of a ferromagnet above which the magnetic moment 
of a material vanishes. Below the Curie temperature T c we can observe 
only the ferromagnetic phase which continuosly disappears at T c . Above T c 
we can observe only the paramagnetic phase. Some phase transitions are 
not in accord with the naively applied classification rules. Thus, the Bose- 
Einstein condensation in an ideal Bose gas is accompanied by a kink in the 
temperature dependence of specific heat at the Bose-Einstein condensation 
temperature but is viewed as a first-order phase transitionPln what follows 
our focus will be on the continuous phase transitions which were studied 
very intensively in the last century, especially in the last four decades. 

Statistical mechanics provides some microscopic models of the continu- 
al 

ous phase transitions. One of such models was invented by E. Ising about 
80 years ago. It became especially well-known after L. OnsageJ^had found 
an exact solution of the model in two dimensions. Let me introduce the Ising 
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model and fix the notations. The model consists of magnetic moments or 
spins which may have two values — \ and \ and which interact with the 
nearest neighbors. The Hamiltonian of the model on a square lattice of 
N x N y — N sites can be written in the form 

Ny 

ix = l i y = l 

where sf x j is the spin variable attached to the site i x ,i y , and Jh and J v 
are the exchange interactions in horizontal and vertical directions, respec- 
tively. If the exchange interaction in Q is positive, the same value of spin 
variables at neighboring sites is favorable, i.e. the exchange interaction is 
ferromagnetic. The spin variable s z may be presented by a half of the Pauli 
matrix 

1 
-1 

and the canonical partition function which determines the thermodynamics 
of the model reads 

Z = ex P i~P H ) = Trex P (~P H ) i ( 2 ) 

where f3 = pp is the inverse temperature. We are also interested in the spin 
correlation functions (sf »• sf a ), where the canonical average means 

x t x ,ly Jx Tjy ' u 

((...)) = |Tr (exp (-j3H) (...))■ 

The spin correlation functions in the limit of infinitely large intersitc dis- 
tances yield the magnetization per site 

which plays the role of the order parameter. Due to the seminal study 
of L. Onsager we understand in great detail the properties of the two- 
dimensional Ising model. 

At zero temperature T = all spins have the same value, say, a, and 
rn z = h . As the temperature becomes nonzero some of spin variables due to 
temperature fluctuations have the value — \ , and thus m z becomes smaller 
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than |. Quantitatively the temperature fluctuations of magnetic moment 
can be characterized by 

2 




At the critical temperature T c the temperature fluctuations completely de- 
stroy the order: the average numbers of spins having the values \ and — | 
are the same and m z — 0. Above T c the average numbers of spins having 
the values | and — ^ remain the same and m z = 0. Analytical calculations 
predict the logarithmic singularity of the specific heat in the vicinity of T c 

c~ln|T-T c |, (4) 

and the vanishing of the order parameter m z while T c — T — > +0 as 

m z ~{T c -T)K (5) 

The correlation length £ which characterizes the long-distance behavior of 
spin correlations, 



exp — , |i — j| — »• oo, (6) 



diverges in the vicinity of T c as 

e-ir-Te]- 1 . (7) 

In one dimension T c = 0. At any infinitesimally small temperature the 
temperature fluctuations destroy the long-range order which exists only at 
zero temperature. The exact solutionSI for the two-point spin correlation 
functions reads 

( s i s l+»> = \ ( tanh = \ 6XP ("| 

i=lncoth^, (8) 

and hence the correlation length £ diverges while T — > T c = as 

C-exp^. (9) 

We do not know the exact solution of the Ising model in three dimensions 
although qualitatevilely the briefly sketched picture for a two-dimensional 
case remains valid. 
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Numerous experimental studies of different substances in the vicinity 
of the continuous phase transition points show that the critical behavior 
is characterized by a set of exponents which may be identical for different 
substances. This remarkable result of universality urged the researchers to 
proceed elaborating scaling concepts, establishing scaling relations and de- 
veloping re-normalization ideas and many systematic renormalization-group 
schemes in order to calculate critical exponents. These issues constitute a 
basic course in t he th eory of classical (i.e. temperature-driven) continuous 
phase transitions!^! 

In what follows I shall not speak about temperature-driven continuous 
phase transitions. I wish to focus on the quantum continuous phase tran- 
sitions. Let me start from a brief discussion of the experiment performed 
by D. Bitko, T. F. Rosenbaum and G. Aeppl0 which demonstrates how 
a phase transition may be driven by entirely quantum rather than tem- 
perature fluctuations. D. Bitko et al carried out their measurements for a 
model magnet lithium holmium fluoride LiHoF4 in the external field H f , 
which turns out to be the experimental realization of the Ising magnet in 
a transverse magnetic field. At low temperatures (below 2 K), the mag- 
netic properties arise owing to the magnetic dipolar interaction of the spins 
of neighboring holmium ions Ho 3+ . These spins prefer to be directed ei- 
ther up or down with respect to a certain crystalline axis and present the 
three-dimensional ferromagnetic Ising model. D. Bitko et al investigated 
the behavior of LiHoF4 as a function of temperature T and an external 
magnetic field if* applied normally to the Ising axis. For this purpose they 
measured the real and the imaginary parts of the magnetic susceptibility 
along the Ising axis at different temperatures establishing T c or if*. Their 
findings are summarized in Fig. ^ 

Let us discuss the experimentally measured phase diagram (Fig. ^| of 
the system which is described by the Hamiltonian of the Ising magnet in a 
transverse magnetic field 

ff = £jeaf^-nX>f, (io) 

hj i 

where s a , a = x,y, z are halves of the Pauli matrices 

-(:;)• -(!o). -G-.). <"> 

JijS are the Ising exchange couplings and f2 is a transverse field. Con- 
sider the behavior of the system as temperature increases at H f = (the 
horizontal axis in Fig. Q). At T = 0, all spins are pointed, say, up (i.e., 
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Fig. 1. The temperature — transverse field phase diagram of LiHoF4. The temperature 
T (the horizontal axis) is measured in K, the (transverse) magnetic field H l (the vertical 
axis) is measured in kOe. Experimental data for the ferromagnetic transition (filled 
circles) are obtained via magnetic susceptibility measurements. 



fully polarized ferromagnetic state). As the temperature becomes nonzero, 
the temperature fluctuations cause some spins of the ground-state config- 
uration to flip, i.e. some spins become pointed down. As the temperature 
increases, the number of flipped spins increases and, at T c = f.53 K, the 
numbers of up and down spins become equal. The numbers of spin up and 
spin down are the same for all temperatures above T c . This is a scenario 
of the conventional temperature-driven continuous phase transition from 
the ferromagnetic to the paramagnetic phase with the Curie temperature 
T c = f.53 K. 

Consider further what happens while we are moving along the vertical 
axis in Fig. Q i-e. while we increase the transverse field H* at zero tem- 
perature T — 0. We immediately observe that the existing order at zero 
temperature may be destroyed in a completely different manner even at 
T = 0. Applying transverse field H we allow tunneling between spin-down 
and spin-up states. And if > H*, the ground state becomes paramag- 
netic even at T = 0. Hence, the quantum fluctuations destroy the order. 
In other words, the phase transition between ferromagnetic and param- 
agnetic phases described by the order parameter -^-X^( s f) 1S driven en- 
tirely by quantum fluctuations. For both T and H l nonzero, D. Bitko et 
al found experimentally a line of continuous phase transitions which sepa- 
rates the ferromagnetic phase (lower region in Fig.^) from the paramagnetic 
phase (upper region in Fig. HJ. For example, at T = O.fOO K they found 
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Hi = 49.3 kOe. Thus, the classical phase transition at if* = and the 
quantum phase transition at T — are connected by a line of continuous 
phase transitions. To develop a correct theoretical description of the equally 
important quantum and temperature fluctuations in this region has been a 
focus of many recent studies. 

For further discussions on this issue as well as for other examples of 
quantum phase transitions see the review article for a general science 
audiencJSTand the bool0of S. Sachdev. 

2. Spin-i Ising Chain in a Transverse Field as the Simplest 
Model for the Quantum Phase Transition Theory 

In early sixties of the last century E. Lieb, T. Schultz and D. MattiJi^ (see 
also the paper by S. Katsure^^) suggested a new exactly solvable model, the 
so-called, spin-i XY chain. In particular case it transforms into the one- 
dimensional spin-i Ising model in a transverse field; this case was examined 
in detail several years later by P. Pfeutyff^ The long-known results for the 
Ising chain in a transverse field may be viewed in the context of the quan- 
tum phase transition theory. The transverse Ising chain is apparently the 
simplest model exhibiting the continuous quantum phase transition which 
can be studied in much detail since many statistical mechanics quantities 
for that model are amenable for rigorous calculations. 

I begin to discuss these results introducing the model. It consists of 
N — > oo spins \ (which are represented by halves of the Pauli matrices) 
which are arranged in a row. Only the neighboring spins interact via the 
Ising exchange interaction. Moreover, the spins interact with an external 
transverse field. The Hamiltonian of the model may be written as follows: 

N N 

ff = 5>< + £j« +1 . (12) 

n— 1 n— 1 

We may impose periodic (cyclic) boundary conditions assuming s^ +1 = s" 
or open (free) boundary conditions assuming the last term in the second 
sum in i|12|) to be zero. The Ising chain without transverse field, fl = 0, 
exhibits the long-range order at zero temperature T = with the order 
parameter (s x ) — X]i( s f) which equals, say, \ (we assume the exchange 
interaction in 112|l to be ferromagnetic, J < 0). The long-range order is 
immediately destroyed, i.e. (s x ) — 0, for any nonzero temperature T > 
0. However, the Ising magnetization (s x ) can be destroyed at T = by 
quantum fluctuations, when we switch on in l|12(l the transverse field f2 ^ 0. 
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Small values of reduce the transverse magnetizations (s x ) and after 
exceeds the critical value f2 c = ^ the Ising magnetization becomes zero. 
That is the scenario which was discussed above in connection with the low- 
temperature properties of LiHoF4. The advantage of the introduced model 
(|12|l is a possibility to follow in great detail the quantum phase transition 
tuned by tt. 

Just after introducing the spin-^ XY chains it was recognized that 
there was an intimate connection between the transverse Ising chain and 
the square-lattice Ising model (see, for example, a redirevation of the On- 
sager solution using the Jordan- Wigner fermionizatioiff^) . The relation ship 
between these models was demonstrated explicitly by M. Suzuki.EllSI Con- 
sider the square-lattice Ising model with the strength of horizontal inter- 
actions Jh > and with the strength of vertical interactions J v > and the 
transverse Ising chain i|12|) with the exchange interaction J < and with 
the transverse field f2. The thermodynamic properties of the square-lattice 
Ising model are equivalent to the groun d-state properties of the transverse 
Ising chain under the relation 

2kT l J l 

Moreover, the temperature-driven continuous phase transition in the 
square-lattice Ising model corresponds to the transverse field-driven con- 
tinuous phase transition at T = in the transverse Ising chain. The critical 
temperature T c corresponds to the critical transverse field Q c - The depen- 
dences such as the Helmholtz free energy, entropy, and specific heat against 
temperature T correspond to the dependences such as the ground-state en- 
ergy, transverse magnetization, and static transverse susceptibility against 
transverse field f2. Finally, the correlation length £ corresponds to the in- 
verse energy gap A, £ ~ A -1 . The behavior of some ground-state quantities 



found l>y P. Pfeutjfi^is tvportrd in Table J and Fiu, 2 



Table 1. Towards the correspondence between thermody- 
namic properties of the square-lattice Ising model and the 
ground-state properties of the transverse Ising chain. 



square-lattice Ising model 


transverse Ising chain 


m z ~ (T c -T)i 
c~ln|T-T c | 
£~ \T-Tc\~ 1 


(s x ) ~ (Q c -n)3 
x z ~ ln|fi- n c | 
a ~ \n - n c \ 
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Fig. 2. The ground-state dependences of the transverse magnetization m z (a), the static 
transverse susceptibility \ z (t>), and the longitudinal magnetization m x (c) on the trans- 
verse field Q for the transverse Ising chain 1121 (J = —2). 

Now I wish to explain briefly how can the statistical mechanics cal- 
culations for the transverse Ising chain be carried out without making 
any approximation. In terms of the spin raising and lowering operators 
s n = S n ^ ls n the Hamiltonian of the transverse Ising chain becomes as 
follows: 

+ ( s n s n+l + s n s n+l + S n S n+1 + S n S n+l) • 

n 

Bearing in mind a case of regularly alternating chain we consider a model 
more general than (|12fl assuming that the Hamiltonian parameters are site- 
dependent, f2 — * fi n , J — > J n — 2In- Although Ijl4(l is a bilinear form 
in terms of s operators, further calculations appear to be complicated 
because of the commutation relations which are of the Fermi type at the 
same site 

{s«,s+} = l, {s~,s~} = {s+,s+} = (15) 
and of the Bose type at different sites 

[ s n > 4J = [ s n, s m] = [4> «m] = °i n ^ Vtl. (16) 

We may use the Jordan- Wigner transformation a to introduce the Fermi 
operators according to the following formulas: 

c n = (-2 S f)(-2 S l)...(-2<_ 1 )^, (17) 

a The Jordan- Wigner transformation of the spin operators to spinless fermions is also 
described in Chapter ? (see Eq. (3) of that Chapter). 
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c+ = (-2 S f)(-2 S |)...(-2<_ 1 ) S + (18) 

Really, the introduced operators always satisfy the Fermi commutation re- 
lations 

{c n , c m | = 5 nm , {c n , c m } — |c n , c m | 0. (19) 

Moreover, the Hamiltonian (|14|) in terms of the Fermi operators (|17fl . (|18|) 

remains to be a bilinear form. Namely, 

n ^ ' 

Y ( C n C n+1 + C n c n+l — c « C n+l — c n c n+l J 

n 

= 7) ^ ] ^71 ~t~ ^ ^ f ^nm'iii "i" 7) i^n^nra^rn C n i? nm C m ) j (20) 



2 / ^ f* / ^ I — Ibtlb lib ' ^ 

n n,m ^ 

where 

A n m — ^n^nm "i" ^ -^n^m,n-(-l ~t~ ~<Y^ n ~ l^rn.,n— 1 = Anni (21) 
^nm 2^ n ^ m i Tl +l 7^ — 1 ^T7i , ?i — 1 -B mn . (22) 

For periodic boundary conditions imposed on (|14fl . the transformed Hamil- 
tonian (|20[) should also contain the so-called boundary term which is omit- 
ted since we send N to infinity. 

To diagonalize a bilinear in Fermi operators form 11201) we perform the 
linear canonical transformation with real coefficients 

Vk (dkn.Cn + hknCn) , 

n 

Vk = 51 (-9 fc « C n + hknCn) ■ (23) 
n 

The resulting Hamiltonian becomes as follows: 

JV 



fc=i 



{%' , Vk" } = <W' {%' - Vk»} = {ny, Vk" } = ( 24 ) 
if the coefficients g^ ni h^n or, more precisely, their linear combinations 

&kn = 9kn + hkn, 
^kn = gkn - hkn (25) 
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11 



satisfy the set of equations 




5> fcj ((A-B)(A + B)) r 




J2*kj((A + B) (A-B)) r 



(26) 



3 



Eqs. iPEt explicitly read 



n n i„-i$ k , n -i + (il + nl - A 2 k ) * fcn + n n+1 i n ^ ktn+1 = o, (27) 



with periodic or open boundary conditions implied. Eqs. (|27|l . I|25|l deter- 
mine the coefficients gkn , hkn in 1|23|) and the elementary excitation energies 
A fc in (33} . 

For the uniform transverse Ising chain, D„ = f2, /„ = /, the bilinear in 
Fermi operators form (|20|l becomes diagonal after performing: i) the Fourier 
transformation 



with k — jj-n and n = — ^, — y + 1, . . . , ^ — 1 if N is even or n — 
— N 2 1 : + 1: ■ ■ • j if ^ is odd and ii) the Bogolyubov transforma- 

tion 




(28) 



+ * i * 4- 



(29) 



with 



i/sin k 



\/ 2A K (A« - e„) ' 




(30) 



and 




= Q + I COSK. 



(31) 
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Knowing the energies of noninteracting fermions Aj, in (|24[1 we imme- 
diately obtain the Helmholtz free energy per site (and hence all thermody- 
namic quantities), 

/ = --L In Trexp (-/?#) 
' ' /?A, 



J_ lnn ^xp^) + exp v 2 



fc=i 

AT 



= In ( 2 cosh 



/3A 



y d/cln^cosh^^ (32) 



(the last line in 1)32(1 refers to the uniform case l|31|) ). Note, that we can ob- 
tain all thermodynamic quantities knowing the distributions of the energies 
of elementary excitations, 

1 N 

fe=i 

dEp{E) = 1, (33) 

or the distribution of the squared energies of elementary excitations 

1 N 



I AE 2 R(E 2 ) = 1. (34) 

j0 



Really, the density of states p(E) (j33) or the density of states R(E 2 ) (f3"4"|l 
immediately yields 



l£d£,(E)ln(2c„ s hf) 
j AEER(E 2 ) In ^2 cosh . 



(35) 



2 

The calculation of the spin correlation functions in the fermionic picture 
looks as follows. First we write down the relations between the spin and 
Fermi operators. We have 

1 1 / J. X / ^ x 1 



c' c 



n - - = ~- (c+ + c„) (c+ - c„) = --<p+<p n (36) 
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where we have introduced the operators ip^ = c+ ± c„ b . Further, 

s n = ^VlPl ■ ■ ■ Pn-lPn-lVn* ( 37 ) 

s l = J^lVl ■■■ <Pn- 1 <Pn- l¥>n • ( 38 ) 



Note that the relations between the operators s*, s}^ and the Fermi op- 
erators are nonlocal since the r.h.s. in Eqs. Ij37(l . I|38|) involve the Fermi 
operators attached to all previous sites. That is in contrast to the similar 
relation (|36J) for the operator s^, the r.h.s. of which contains only the Fermi 
operators attached to the same site n. 

Now we can rewrite the spin correlation functions in fermionic language. 
For example, 

(«+m) = \{<Pt<Pn<Pn+mVn+m)> ( 39 ) 
(«+m) = \{<Pn<P$+l<Pn+l ' ' ' Pt+m-lVn+m-lVn+m) • ( 40 ) 

To get the latter result we use the relations 

{<Pt><Pm} = - {^n.Vm} = 2S nm, {^,9m} = °- ( 41 ) 

The operators ip^ are linear combinations of operators T)k, Tjt involved 
into (EU), 

N 

p=l 
N 

or for the uniform chain 

¥>« = ^= X! exp ( Ik + f«) (^k + V-") ' 



b The operators are related to the Clifford operators T 1 , T 2 (Chapter ?, Eq. (5)) as 
follows: tpX = T^, ipn = -il" 2 . 
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Therefore, we calculate <|39[) . (I4U|) using the Wick-Bloch-de Dominicis the- 
orem . Namely, 



4(«+m) = (<p£<Pn<P£+m<Pn+m) 
= (<Pt<Pn){<Pt+m<Pn+m) 
III ><¥> II > + (ft fn+ mi vrn ml ' 



(44) 



The r.h.s. of Eq. (|44|l may be compactly written as the Pfafnan of the 4x4 
antisymmetric matrix 



A(s z s z , ) 



Pf 



/ ififn) 

-(<Pt<Pn) 

\-{<P£Pn+m) -{<Pn<Pn+m) 



ml Vrn 
(vt+mVn+m) 



(45) 



+7TJ / 



/ 



Similarly, 



= Pf 



/ (<PnV£+l) (<Pn<Pn+l) 



V-^nVifm) -(^n+l^n+m) -(^n+lVn+m) 



o / 



(46) 



The elementary contractions involved into l|45|l . (|46|l read 



N 



(rifl) = J! ^kn^km = 5 nm , 



k=l 



N 



(Pnfm) = ^kn^km tanh 



k=l 
N 



(fn Vm) = - X] ^kn^km tanh 



2 



k=l 
N 



(47) 



fe=i 



c See also Section 3.1.2 in Chapter ?. 
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In the uniform case, instead of (|47|l we have 



N 



cxp (in (n — m) ) 



fi + /cxp(-iK) tanh /3A K 



(<Pn <Pm) = eXp ( n - m )) 



Q + /exp (— i/t) 
AT 



tanh 



2 

/3A„ 



(48) 



It is worthy to recall some properties of the Pfaffians which are used in 
calculating them. In the first numerical studies the authors used the relation 

(PfA) 2 = detA (49) 

and computed numerically the determinants which gave the va lues of Pfaffi- 
ans. On the other hand, the Pfaffian may be computed directly, ^ ^ ^ ^ noting 
that 

Pf(U T AU) =dctUPfA 

and that 

/ R 12 ...0\ 
-i?i2 ... 
R 34 ■ ■ ■ 

-#34 ... 



(50) 



Pf 



12^1.34 



(51) 



\ ...0/ 

Let me finally s how how the results for the transverse Ising chain ob- 
tained by P. Pfeut}E2 

arise in the described approach. Recalling Eq. (|24l) 
(or Eq. 1321) ') it is not hard to see that the ground-state energy per site can 
be written as 



eo 



2vr 



dn - 



A„ 



1 PTT 

= / dny/VL 2 + 2QIcosk + P 

= / d«VA 2 + 2Acosk + 1 

1(1 + X) 



\ 



2%/A 
1 + A 



f (1 + A) ( 7T gVX 
tt 2' 1 + A 



(52) 
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16 O. Derzhko 

where A = y and 

^(|,o) = Jj d<p\J 1 - a 2 sin 2 
is the complete elliptic integral of the second kind !^ In the vicinity of a = 1 

(see Ref. E01 0} and therefore in the vicinity of A c = 1 the ground-state 
energy 152J1 contains the nonanalytic contribution 

e ~ (A-A c ) 2 ln|A- A c | . (54) 

As a result the ground-state transverse magnetization contains the nonan- 
alytic contribution 

m z ~ (A- A c )ln|A- A c | (55) 

and the ground-state static transverse susceptibility exhibits a logarithmic 
singularity 

X 2 ~ln|A-A c |. (56) 

At A c = 1 the energy spectrum is gapless, i.e. there is a zero-energy ele- 
mentary excitation. In the vicinity of A c = 1 the energy gap is given by the 
smallest value of A K = \/£l 2 + cos k + 1 2 and hence 

A~|A-A C |. (57) 

To end up the discussion of the ground-state properties of the transverse 
Ising chain, let me emphasize that the quantum phase transition appears 
since the spin variables are the q- numbers rather than the c- numbers. To 
demonstrate this explicitly we may consider spin-vector (instead of spin- 
matrix) transverse Ising chain. The model consists of N — > oo 3-componct 
vectors 

s = (s x , s v ,s z ) — (s sin cos cj>, s sin 8 sin (j>, s cos 6) 
which are governed by the Hamiltonian 

n n 

Here s is the value of the spin which plays only a quantitative role and 
further is put s = |. To obtain the ground-state energy we should place 
all spins in xz plane putting <j) n — (I < 0). Moreover, the angles n must 
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Quantum Phase Transitions in Alternating Transverse Ising Chains 17 

minimize the sum of contributions due to the interaction with the field and 
due to the intersite interaction. Thus, the ground-state energy ansatz reads 

E (8) = ^Nfl cos 6- ^N\I\ sin 2 6 (59) 

where is determined to minimize Eq{&) l|59|) . Obviously 

( 1, if O < -2|7|, 
cos<9= I -2]7T»if -2|7| <ft<2|7|, (60) 

[ -l, if 2\i\ < n. 

As a result, we find that the ground-state energy per site is given by 

r §,if n<-2|/|, 

e = I -2|/|<n<2|/|, (61) 

{ if 2|7|<ft. 

Moreover, the x- and z-magnetizations which can be obtained after inserting 
(|fJU)l into the formulas m x — ssin9 and m z — scos9, respectively, behave 
as 

f 0, if n < -2|I|, 

m* = l i^i if -2|7|<n<2|/|, (62) 

[ o, if 2\i\ < n 

and 

f |, if n<-2|i|, 
mZ = { -wr [i -2|/|<n<2|j|, (63) 
[ -|, if 2|/|<n, 

respectively. Evidently, there is no singularity in the ground-state depen- 
dence of the static transverse susceptibility \ z = on the transverse field 
fl. As the transverse field increases from zero to 2|7|, the on-site magnetic 
moments smoothly change their direction from x- to z-axis. The quantum 
phase transition appears only due to the quantum fluctuations which arise 
when the spin components do not commute and no phase transition oc- 
curs when the spin components are classical variables (compare Fig. and 
Fig. El. 

3. Spin-i Ising Chain in a Transverse Field with Regularly 
Alternating Hamiltonian Parameters: Continued 
Fraction Approach 

We are turning to a discussion of the effects of regularly alternating ex- 
change interactions and fields on the quantum phase transition inherent in 
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the transverse Ising chain. Note, that there was a great deal of work done 
examining the effects of different modifications of the skeleton model on the 
quantum phase transition. Thus, J. M. LucH^ 1 analyzed the critical behav- 
ior of the chain with an aperiodic sequence of interactions, D. S. FisheJ22 
performed an exte nsiv e real-space renormalization group study of random 
chains, F. Igloi et aP^ reported a renormalization-group study of aperiodic 
chains. However, a simpler case of regularly alternating transverse Ising 
chain still contains a good deal of unexplored physics. It appears possible 
to obtain analytically exact results for thermodynamic quantities of such 
models using the Jordan- Wigner fermionization and the continued fraction 
approach ! ' ^ ' Moreover, for certain sets of the Hamiltonian parameters 
the exact results may be obtained for the spin correlation functions!^ In 
all other cases, spin correlations may be explored numerically at a high 
prec ision l evel considering sufficiently long chains up to a few t housand 
sites™ It is appropriate to mention here the old studied^HEHl referring 
to the chains of period 2 which were exte nded fu rther to the one-dimensional 
anisotropic XY model on superlattices!^^ The quantum critical points 
in the anisotropic XY chains in a transverse field with periodically vary- 
ing intersite interactions (having periods 2 and 3) have been determined 
recently!*^ using the transfer matrix method. 

I begin with explaining how the thermodynamic quantities can be de- 
rived using continued fractions. For this purpose let me recall that the en- 
ergies of the Jordan- Wigner fermions which determine the thermodynamic 
properties of a spin chain are involved into Eqs. 127(1 . Consider, for example, 
the first set of equations in l|27|) which can be rewritten in the matrix form 
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Quantum Phase Transitions in Alternating Transverse Ising Chains 19 



as follows: 



where 



H 



(H-Afl) 



$fe2 
$fe3 



= 



+ 



o 







V ! 



(64) 



(65) 



and 1 denotes the unit matrix. Thus, all A 2 , are the eigenvalues of the 
N x N matrix JBSJ entering To find their distribution R(E 2 ) we 
may use the Green function approach. Consider a matrix G composed of 
the elements G nm — G nm {E 2 ) which is introduced by the equation 



(£ 2 1-H)G = 1. 



(66) 



Having composed a matrix U + of the eigenvectors of the matrix H 16511 . 



UHIT 



/Af ... \ 
' Af . . . 



V ...A% J 



one finds that 



( E l — Af 



£ 2 - A^ 



V o 



\ 





.E 2 -A%J 



UGU+ = 1. 



As a result 



N N 1 

TVG = £(UGU+) fefe =£ * 



fc=i 



(67) 



(68) 



(69) 
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20 O. Derzhko 

Making the substitution 

E 2 -> E 2 ± ie, e -► +0 (70) 
and using the symbolic identity 

E2 _l l± . ie =^^MT^(E^-Al), e^ + (71) 
one arrives at the relation 

N N 
n— 1 n— 1 

Thus, our task is to find the diagonal Green functions G nn defined by Eqs. 

(El- 
Considering the equation for G nn , 

{E 2 — f2 n — I n —i) G nn — Q n -iIn-\G n -i^ n — Q n I n G n +i jri = 1, (73) 

and then the equations 

(E 2 — ^ n _i _ In-2) G>i-l,fi — f^n— lln— 2G n -2,n — ^n-l^n-lG Im = 0, 
(E 2 — &n-2 ~ In-3) Gn-2,n ~ ^n-3ln-sG n -3 >n ~ ^n-2ln-2G n -l, n = 

(74) 

etc. to determine G ^ 1 - n | etc., and the equations 

(E 2 — n n+1 — J n ) G„ + i „ — Q n I n G nn — Q n +lIn+lG n +2,n = 0, 
(_B 2 — f2 n+ 2 — I n +l) G n +2,n — fin4-liiri+lCn+l,Ti — ri n +2^ri+2G rl -|-3 :Tl = 

(75) 

G G 

etc. to determine , G "*°'" etc., we easily find the following continued 

fraction representation for the diagonal Green functions 

G 1 



a: 



E 2 -ni-i 2 _ 1 -A 

O 2 T 2 



2 _ A - _ A+ ' 



02 12 

P2 _ O 2 - f 2 "™-2 J „-2 

^ "n-l i n-2 E 2 -n 2 



ft 2 I 2 

A+ = n , J3 ■ (76) 
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A crucial simplification occurs if a sequence of the Hamiltonian param- 
eters is periodic, 

Then the continued fractions in (|76|l become periodic and can be evaluated 
exactly by solving a square equation. Consider, for example, the case p = 1, 
i.e. the uniform chain. Then A~ = A~ = A (do not confuse with the energy 
gap A) and A satisfies the following equation 

n 2 i 2 

A - EP-w-r-A - (77) 

The solution of Eq. (|77|) . 



A = - [E 2 - fl 2 - I 2 ± y/(E 2 -Q 2 - I 2 ) 2 - 4f2 2 / 2 J , (78) 



G nn = T , (79) 

\2 



yields 



'{E 2 -tt 2 - i 2 y -AWl 2 
and hence in accordance with l|72|) 

[—> • y,if (n - if < e 2 < (n + if , 

R(E 2 ) = I V 4 "^ 2 -^ 2 -" 2 -^) 2 (80) 
[ 0, otherwise. 

We may easily repeat the calculations for p = 2 or p = 3. The desired 
density of states in these cases is given by 

( \z p ^(e 2 )\ - f A (W 2\ ^ n 

^)J pV^^ lf ^)>o. (81) 

( 0, otherwise, 

where Z p _i{E 2 ) and *42p(-E- 2 ) = — TifLi (E 2 — a j) are polynomials of the 
order p — 1 and 2p, respectively, and aj > are the roots of A2 P (E 2 ). 
Explicitly, 

Zi(£ 2 ) = 2E 2 -nj-n 2 ,- l\ - 1$, 

A A {E 2 ) = 40?0|J?J| 

- (E 4 - (nl + + if + if) e 2 + njnl + i 2 ijf 

= -(E 2 - ai ){E 2 - a 2 ){E 2 - a 3 )(E 2 - a 4 ), 

{a J } = ^(n 2 1 + n 2 2 + i 2 + i 2 
±J(n 2 + n 2 + 1 2 + 1 2 ) 2 - MShSh ± hh) 2 ] } ; (82) 
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22 O. Derzhko 



Z 2 (E 2 ) = 3E 4 - 2 (pi + tt\ + fl 2 3 + I 2 + I 2 + I 2 ) E 2 

+ nfif + zfif + njnj + npi + ip 2 
+n 2 n 2 + n 2 i 2 + i 2 i 2 1 

- (e 6 - (n? + nl + n§ + 2? + / 2 2 + if) e* 
+ (njnl + + ijfif + njnl + n|jf + jf jf 
+n 2 n 2 + rtpl + i 2 i 2 ) e 2 



2 ^ ^2 ^3 



i 1 , 

= -(E 2 - ai )(E 2 - a 2 )(E 2 - a 3 )(E 2 - a 4 )(E 2 - a 5 )(E 2 - a 6 ), (83) 

where ai, . . . , ag are the solutions of two qubic equations which follow from 
the equation Aq(E 2 ) = 0. 

Knowing the density of states R(E 2 ) l|81|l we immediately obtain all 
quantities of interest. Thus, the gap in the energy spectrum of the spin 
model is given by the square root of the smallest root of the polynomial 
A2 P (E 2 ). Further, the ground-state energy per site is given by 



e = - / dEE 2 R(E 2 ). (84) 
Ja 

Recalling Eq. (|35|l for the Helmholtz free energy per site we find that the 
specific heat per site is given by 

r°° ( M \ 2 



* = 2 /„ iEER{E \^£m) ■ (85 » 

We assume that Sl n = ft + Afl„ and define the transverse magnetization 
per site and the static transverse susceptibility per site as follows: 

= % < 86 > 



and 

dm z 



(87) 



Unfortunately, the elaborated approach yields only the density of states 
(|34|l and thus is restricted to the thermodynamic quantities. To obtain 
the spin correlation functions (s"s" +n ), a = x, z of the regularly al- 
ternating transverse Ising chain l|4^|l . l|4Tjl we use the numerical 
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approach.™™ The on-site magnetizations can be calculated from 
the spin correlation functions using the relation 

m>^lim_( S ; S « +rp ), (88) 

where 1 <C ji <C N is one of the p consecutive numbers taken sufficiently 
far from the ends of the chain and ji— j\ = 0, 1, . . . ,p — 1. Assuming that 

(.?^>-(.?x.w~«p(-p) m 

for large r (1 <C j, j + rp <C N) we can determine the correlation length 

4. Effects of Regularly Alternating Bonds/Fields on the 
Quantum Phase Transition 

Now we shall use the analytical results for the ground-state and thermody- 
namic quantities and the numerical data for the spin correlation functions 
to discuss the effects of regularly alternating Hamiltonian parameters on 
the quantum phase transition inherent in the transverse Ising chain. 

We start with the energy gap A. The quantum phase transition point is 
determined by the condition A = 0. The density of states R(E 2 ) permits us 
to find the energy gap A since, as it was mentioned before, the smallest root 
of the polynomial A2 P {E 2 ) (see is the smallest elementary excitation 

energy squared. Therefore, the energy spectrum becomes gapless when 

A 2p (0) = 0. (90) 

In the case of period 2 the condition |Jj5DJ, A±{Q) = 0, yields 

n 1 f2 2 = ±/i/ 2 . (91) 

In the case of period 3 the condition (j^DJ, Ae(0) = 0, yields 

Q 1 Q 2 il 3 = ±11/2/3- (92) 

It should be emphasized here that we have rederived the long- 
known condition of the zero-energy elementary excitations obtained by 
P. Pfeuty!^ P. Pfeuty showed that a nonuniform transverse Ising chain 
becomes gapless if 

tt 1 n 2 S} 3 ...tt N =±hhh ...In- (93) 

(In fact, Eq. (6) of the P. Pfeuty's papeJ^does not contain two signs; minus 
immediately follows from symmetry arguments (after performing simple 
rotations of spin axes) and is important for what follows.) 
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24 O. Derzhko 

Analysing the conditions for quantum phase transition points (|91|l , (|92|) 
we immediately observe that a number of quantum phase transition points 
for a given period of nonuniformity is strongly conditioned by specific values 
of the Hamiltonian parameters. Consider, for example, the case p = 2 and 
assume fii 2 = fi ± Ail, Ail > 0. Then, for a small strength of transverse- 
field nonuniformity Ail < y/\ I1-Z2] the system exhibits two quantum phase 
transition points il c — ±-y/Af2 2 + |Jii2|- Moreover, the ferromagnetic phase 
occurs for \il\ < \J Ail 2 + \I\l2\ whereas the paramagnetic phase occurs 
for J Ail 2 + I/1/2I < |^|- If a strength of transverse- field nonuniformity is 
large enough, Ail > ^/|ii/ 2 |, the system exhibits four quantum phase tran- 
sition points il c = Ail 2 ± |/i/2|- Moreover, the ferromagnetic phase oc- 
curs for \J Ail 2 — I-Z1/2I < |f2| < \J Ail 2 + I/1/2I whereas the paramagnetic 
phase occurs for low fields, \il\ < J Ail 2 — I-/1/2I, and for strong fields, 
J Ail 2 + I/1/2I < \Q\- In the case Ail = ^ii^l we get three quantum 
phase transition points, il c = ± a/2 | ^1^2!, . The fields ±\/2|A^2| corre- 
spond to the transition between ferromagnetic and paramagnetic phases. 
At il = only a weak singularity occurs. We shall motivate the statements 
about the phases which occur as il varies considering the dependence of the 
ground-state Ising magnetization on il (see below) . The important message 
which follows from this simple analysis is that the number of quantum phase 
transitions in the regularly alternating transverse Ising chain of a certain 
period strongly depends on a specific set of the Hamiltonian parameters. 

Let us note that in our treatment we assume il n = il + Ail n , fix Ail n 
and /„, and consider the changes in the ground-state properties as il varies, 
thus, breaking a symmetry between the transverse fields and the exchange 
interactions. Alternatively, we may assume I n = I + AI n , fix il n and AI n , 
and assume / to be a free parameter. In this case the quantum phase tran- 
sition is tuned by varying /. Naturally, in general, the quantum phase tran- 
sition conditions l|91l) . (I92|) may be tuned by some parameter (s) effecting 
the on-site fields and the intersite interactions. 

Let us consider a critical behavior of the system in question. To get the 
behavior of the ground-state properties in the vicinity of the critical fields 
we start with expanding the smallest root of the polynomial A2 P (E 2 ) with 
respect to deviation of the field from its critical value e = il — il c . For p = 2 
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with fli 2 = fi ± A17 and il c = Aft 2 ± \I\I2\ we have 

ai = il 2 + Ail 2 + ~ (I 2 + II) 



(ft 2 - ft 



2\ 



2 (ft 2 + Aft 2 + i (i? + if)) ' 
As a result, the energy gap A = ^/ai decays linearly 



(94) 



if ft c and as 



A~|c| (95) 
A ~ e 2 (96) 



if ft c = 0. The latter case occurs when Ail = y/Jlvhl- 

Now we can evaluate the ground-state energy 1(84(1 in the vicinity of 
the critical point. The nonanalytic contribution in the r.h.s. of Eq. 1(84(1 is 
coming from the energy interval between y/ai (which is proportional either 
to |e| or to e 2 ) and y/cia (a± < a 2 < a 3 <...); all the other intervals of 
energies yield only analytical contributions to the ground-state energy eo 
with respect to e. Therefore, we have 

^ r\fa~2 E 2 f(E 2 ) 

eo = / H E — + analytical with respect to e 2 terms 

27r J^ai VE 2 - ai 

~ e 2 In |e| + analytical with respect to e 2 terms. 

(97) 

(Here f(E 2 ), f(0)^0 is some function the explicit expression of which is 
not important for deriving the asymptotic behavior.) Eq. 1(97(1 is valid for 
ft c 7^ when ^/al ~ |e|. As a result 

m z ~ eln |e| + analytical with respect to e 2 terms, (98) 

and 

X z ~ In |e| + analytical with respect to e 2 terms. (99) 

Thus, the critical behavior remains unchanged in comparison with that of 
the uniform chain (compare with (|54Jl . 1(55(1 . 156(1 . ((57(1 ). This can be also 
nicely seen in Fig. 0] where some results referring to the chains of period 2 
with Ii = I2 = 1 and Ail = 0.5 and Ail = 1.5 are collected. 
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-2.0 0.0 2.0 n -2.0 0.0 2.0 n 



Fig. 4. Towards the ground-state properties of the transverse Ising chain of period 2, 

i x = i 2 = i, n i>2 = n ± An, An = 0.5 (n c = ±VT^>) (a - c ), An = 1.5 (n c = 

zbVl-25, ±\/3.25) (d - f). The dependences of energy gap (dotted curves in (a), (d)), 
transverse magnetization m z (solid curves in (a), (d)), static transverse susceptibility 
\ z (dashed curves in (a), (d)), longitudinal sublatticc magnetizations |mj| (solid curves 
with up and down triangles in (b), (e) which are obtained from the data for N = 600), 
and correlation length £ x (triangles and squares in (c), (f) correspond to the data for 
N = 300 and N = 600, respectively), on the transverse field Q. 
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Let us turn now to the case O c = which occurs for the chain of 
period 2 with AQ = ^/\I\h\ - In this case the smallest root a\ tends to zero 
proportionally to e 4 as Q approaches £l c and as a result 

e ~ e 4 In |e| + analytical terms (100) 

and, therefore, 

m z <~ e 3 In |e| + analytical terms (101) 

and 

X z ~ e 2 In |e| + analytical terms. (102) 

Thus, the transverse susceptibility y z is finite in the critical point Q c = 
and only its second derivative with respect to the field exhibits a 

logarithmic singularity. This weak singularity may be called the fourth- 
order quantum phase transition in the Ehrcnfest sense. 

I proceed with exact analytical results for a regularly alternating trans- 
verse Ising chain showing how the ground-state wave function |GS) can 
be obtained exactly for special sets of the Hamiltonian parameters. Let us 
consider the chain of period 2 with 0^2 = f2 ± Af2. It is obvious that 

|GS) = ...| T), | T), + i.. ■ (103) 

as f2 — > — oo and 

|GS) = ...||),||) j+1 ... (104) 
as f2 oo. We can also expect that at O = 

|GS) = ...||),|T), + i... (105) 

if Aft > |/i|,|/ 2 | and 

ICS) = . . . (| [) j ± | ±= (| ± | T>i+i) • • • (106) 

if Ail = 0. Less evident is the ground-state wave function |GS) for £1 = 
=FAfi when fti = 0, Q 2 = -2Aft or fti = 2Aft, ft 2 = 0. To obtain the 
ground-state wave functions for these values of the transverse field let us 
note that after performing the unitary transformations 

U = ... exp (ms x n s v n+l ) cxp (}TTs x n+1 s v n+2 ) . . . (107) 

and 

^ = ...cxp(i|<)cxp(i|< +1 )... (108) 
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the Hamiltonian of the transverse Ising chains with the fields f2 1 f2 2 ■ ■ ■ an d 
the interactions ■ ■ ■ transforms (with the accuracy to the boundary 
terms) into the Hamiltonian of the transverse Ising chains with the fields 
■ ■ ■ and the interactions ftift 2 ■ ■ . , 



R Z UHU + R Z+ =J2 T i 



(109) 



For ft = =fA51 the Hamiltonian R Z UHU + R Z represents a system of nonin- 
teracting clusters which consist of two sites. We can easily find the ground- 
state of a two-site cluster with the Hamiltonian (ft = — Aft) 



#12 =hsf +I 2 s|-4Afis?s§. 



(110) 



It reads (hh > 0) 



Cl 



C2 



|GS) 12 =ciU)iU) 2 
1 I + V^ 



•Pa|t>i|T>s 



Aft 2 



V2 v 7 / 2 + UP + Aft 2 + Aft 2 ' 

1 Aft 
V2 a// 2 + IV/ 2 + Aft 2 + Aft 2 ' 

I=\{h+h)- 



(111) 



As a result, the desired ground-state wave function reads 

|GS) = U+R z+ . . . (d| |) n+1 | |>„ +2 + pa] T) n +i| T>«+ 2 ) 

• (Cl| |)n+3| l)n+4 + C 2 | T)n+ 3 | T)n+ 4 ) .... (112) 

For ft = Aft, the ground-state is again given by (|112fl . (|lll|l with the 
change n — > n — 1, Aft — > —Aft. 

Knowing the ground-state wave function we can easily calculate different 
spin correlation functions. For example, for ft = —Aft according to Eq. 
(ITT2ft we get 



^n+l°n+2l 



\ S n+l s n+3/ 



3 n+l c 'n+4/ 



1 



(ci " C?) 



'n+2' : 'ri+6/ 



(113) 
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( s n+l s ri+2) = ( s n+l s n+4) = ■ ■ ■ 
= ( s ri+2 s n+3) = ( s n+2 s n+5) = • • • 
— { s n+2 s n+i) = { s n+2 s n+&) = ■ • • = 0, 

«+i< +3 ) = «+i<+5> = • • • = (cic 2 ) 2 • (H4) 
For Q = Af2, the correlation functions follow from (|113fl . (|114|l after the 
change n-tn-l. 

Now we discuss the exact numerical results for finite chains presented 
in Figs. 0] 03 The ground-state Ising (longitudinal) sublattice magnetiza- 




-0.2 o.o 0.2 n 



Fig. 5. Towards the ground-state properties of the transverse Ising chain of period 2, 
I 1 = I 2 = l, Sly = n± AQ, AO. = 0.99 (open symbols), AQ = 1.01 (full symbols). The 
longitudinal sublattice magnetization |mj| (a) and the correlation length £ x (b) against 
the transverse field Q for chains having 300, 600 and 900 sites (triangles, squares and 
circles). 

tions |mj| indicate different phases and phase transitions between them. 
Their dependence on the transverse field are shown in Figs. 01 for 
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Aft < y/\hh} (Fig. Hd and Fig. EJi (open symbols)) and An > y/\hh\ 
(Fig. 0k and Fig.[SjL (full symbols)). In accordance with the analytical re- 
sults for A, m z and \ z the numerical data show the existence of either 
two phases (quantum Ising (ferromagnetic) phase for < -J Ail 2 + \I\l2\ 
and strong-field quantum paramagnetic phase otherwise) or three phases 
(low-field quantum paramagnetic phase for |f2| < y/ AO 2 — 1 7i X2 ] , quantum 
Ising (ferromagnetic) phase for -J AO 2 — I-Z1/2I < |0| < -J Aft 2 + I/1/2I and 
strong-field quantum paramagnetic phase otherwise). The longitudinal on- 
site magnetizations mj are nonzero in quantum Ising phases and become 
strictly zero in quantum paramagnetic phases. The transverse magnetiza- 
tion m z in quantum paramagnetic phases is almost constant being in the 
vicinity of zero in the low-field phase and of saturation value in the strong- 
field phase, thus producing plateau-like steps in the dependence m z versus 
Q (see the solid curve in Fig.0Ji). The correlation length f x (Figs-EJijEfjEb 
(full symbols)) illustrates that the transition between different phases is ac- 
companied by the divergency of £ x . Comparing the data for N — 600 and 
N = 900 reported in Fig. 03 one observes a strong size-dependence of the 
correlation length value £ x about the critical fields (full squares and circles 
in Fig. 03 ) and a weak size-dependence of £ x for other values of fl (for 
example, open squares and circles in Fig. 0d) . It is interesting to note that 
short chains (N — 20) are already sufficient to reproduce a correct order- 
parameter behavior in the quantum Ising phase away from the quantum 
critical point, but not in the quantum paramagnetic phases!^ The chain 
length of up to N — 900 sites is clearly sufficient to see a sharp transition 
in the order parameter (Figs.0jD,0J;) and even to extract reliable results for 
£ x from the long-distance behavior of (sjsj +n ) as can be seen from the data 
presented in Figs. 0Ji, 0f and in Fig. [SJx We also note that the numerical 
data for |mj| and £ K at n — =fA57 coincide with the exact expression (II 1311 . 
In particular, the values of the longitudinal sublattice magnetizations for 
these fields are ^ and \ \c^ — c\\. 

Finally, the low-temperature dependence of the specific heat c at dif- 
ferent fl (> 0) obtained from the exact analytical expression (|85fl confirms 
the existence of either two phase transitions (Fig. Et) or f° ur phase transi- 
tions (Fig.^) depending on the relationship between Att and y^/ii^l- At 
the quantum phase transition points the spin chain becomes gapless and 
c depends linearly on T. Really, at the quantum phase transition point as 
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Fig. 6. The low-temperature dependence of the specific heat for the transverse Ising 
chain of period 2, h = J 2 = 1, fti,2 = fi ± AC, AC = 0.5 (a), AC = 1.5 (b). 

T ^ 

^^r d ^(2i) 2+2 r™ 2) (^) 2 - 

(115) 

(Here C is some constant the value of which is not important for deriving 
the asymptotic behavior.) The second term in (|115f) disappears in the limit 
T (E> kT) and as a result 

c~T. (116) 

The ridges seen in Figs.[Sli,[^j3 correspond to the maxima in the depen- 
dence c versus fi as T varies. They single out the boundaries of quantum 
critical regionsP These boundaries correspond to a relation A ~ kT that 
can be checked by comparison with the data for A versus Q reported in 
Figs. EH, Eli- As can be seen from Fig. the c(T) behavior for £1 slightly 
above or below f2 c changes crossing the boundaries of quantum critical re- 
gions. Furthermore, we notice that for AJ7 = 1.5 (Fig. Eh) an additional 
low-temperature peak appears in the temperature dependence of the spe- 
cific heat. 

To end up, let us turn to a chain of period 3 for which the analytical and 
the numerical calculations presented above can be repeated. The quantum 
phase transition points follow from the condition H92J) . Depending on the 
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parameters, either two, four, or six quantum phase transitions are possible. 
This behavior is illustrated in Fig. for a chain of period 3 with h = h = 



Fig. 7. Towards the number of quantum phase transitions present in the transverse 
Ising chain of period 3, h = h = h = b ^1,2,3 = H + AS7i >2 ,3, Af2i + Af2 2 + AS7 3 = 0. 
The dark, gray, or light regions correspond to parameters for which the system exhibits 
two, four, or six quantum phase transitions, respectively. At the boundaries between 
different regions weak singularities occur. 

J 3 = 1, fi n = O + Afi„, AO x + AS1 2 + AS1 3 = 0. Moreover, such chains 
may exhibit weak singularities. For example, for the set of parameters at 
the boundary between dark and gray regions there are three critical fields; 
at one of them a weak singularity occurs. 

To summarize, I have shown how the exact results for the thermodynam- 
ics of a regularly alternating spin-i Ising chain in a transverse field can be 
derived using the Jordan- Wigner transformation and continued fractions. 
Furthermore, for special parameter values the exact ground-state wave func- 
tion and spin correlation functions can be obtained as well. For parameters 
for which the correlation functions are not accessible for rigorous analytical 
analysis they can be obtained from exact numerical results for finite but 
large systems. From these numerical data we can extract the order param- 
eter and the correlation length in high precision. We have found, that the 
quantum phase transition occurring in the uniform transverse Ising chain 
is not suppressed by deviation from uniformity in the form of a regular al- 
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ternation of bonds and fields. On the contrary, field alternation may lead to 
the appearance of additional quantum phase transitions tuned by the field. 
The number of phase transitions for a given period of alternation strongly 
depends on the precise values of the parameters of the model. We have ex- 
amined the critical behavior of the energy gap A, the ground- state energy 
eo, the transverse magnetization to 2 , and the static transverse susceptibility 
X z and have found the same critical indices both for a nonuniform and uni- 
form chain. For some sets of the Hamiltonian parameters, the ground-state 
quantities may exhibit a weak singularity. 



5. Related Models 

In the remainder of this chapter I shall discuss another simple models which, 
apart from the transverse Ising chain, exhibit a quantum phase transition. 
Namely, I consider the spin-i isotropic XY (i.e. XX) chain in a transverse 
field and the spin-i XYZ chain in an external magnetic field directed 
along z axis focusing on the effects of regularly alternating Hamiltonian 
parameters on the dependence of the ground-state quantities on the external 
field. 

The uniform transverse XX chain described by the Hamiltonian 

N N 
n— 1 n—1 

by employing the Jordan- Wigner transformation l|17fl , (|18fl can be mapped 
onto the system of noninteracting spinless fermions with the Hamiltonian 

tf = 5>«(c+c„-i) (118) 

where the elementary excitation energies A K are given by 

A k = Q + 2Icosk. (119) 
Alternatively, we may find the density of states p(E) (|33() which is given by 

f— — 1 =, if ai 2 - (e - nf > 0, 

p (E) = I ?V« 2 -(£-n) 2 ' (120) 
[ 0, otherwise. 

As a result, the thermodynamic properties of the spin chain 1117(1 can be 
easily analyzed. 

In contrast to the transverse Ising chain (and to any chain with 
anisotropic XY interaction in a transverse field) Yln=i s « commutes with 
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the Hamiltonian describing the isotropic XY interspin interaction and, 
therefore, a change in the transverse field SI has a different effect on the 
ground-state quantities of the spin chain (|117|l . As can be seen from (|118Jl . 
(111911 the transverse field plays a role of the chemical potential controlling 
a filling of the fermion band. The transverse XX chain remains gapless 
until |S1| < 2\I\. If |S1| exceeds 2\I\, the energy gap A opens linearly. This 
produces singularities in the ground-state quantities. For example, a square- 
root singularity in the zero-temperature dependence of the static transverse 
susceptibility x z on the transverse field SI. 

After introducing a regular alternation, i.e. after making a substitution 
in l|117|) SI — ► Sl„, / — > I n with a periodic sequence of parameters 

O1/1SI2/2 • ■ ■ Slp/pSli/iSl2^2 • ■ ■ Qplp ■ ■ ■ , 

we can easily obtain with the help of continued fractions the density of 
states p{E) of the Jordan- Wigner fermions which represent the regularly 
alternating transverse XX chain, and hence analyze the effects of regularly 
alternating Hamiltonian para meters on the thermodynamic properties of 
the considered spin systenfSl (for another approach see Ref . 35 ) . The main 
consequence of the introduced periodic nonuniformity is a splitting of the 
initial fermion band into several subbands the number of which does not 
exceed the period of the chain p (for special (symmetric) values of the 
Hamiltonian parameters one may observe a smaller than p number of sub- 
bands). The transverse field SI (Sl„ = CI + AS1„) controls a filling of the 
fermion subbands. Again the energy gap A disappears/appears linearly 
and the critical behavior remains as for the uniform chain. Note, that the 
ground-state dependence of the transverse magnetization 

1 r°° E 
m z = --J dEp(E)t<mh — (121) 

on S7 is composed of sharply increasing parts separated by horizontal parts 
(plateaus) in accordance with a famous conjecture of M. Oshikawa, M. Ya- 
manaka and I. Affleck (see Ref . . 

A more general situation is the case of regularly alternating XYZ chain 
in a magnetic field. Let us consider the Hamiltonian 

H = IK + E (»n +1 + + I Z « + l) , 

n n 

I* = l+j+(-l) n 5, iv = l- j 5, (122) 

which captures the interplay between an exchange interaction anisotropy 
(7, I z ), exchange interaction modulation or, more precisely, exchange in- 
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teraction dimerization (5), and the uniform magnetic field (Q). Such a 
model has been studied recently in some details- using the bosonization 
and the Lanczos diagonalization technique. Let us first 
discuss the free fermion case I z — 0. After the Jordan- Wigner transforma- 
tion (|17(1 . (|18(l the Hamiltonian can be readily diagonalized (see Refs. GBl 
1291 1411 l3*?j) . The critical field values at which the spin system 1)122(1 becomes 
gapless are given by 

n c = ±^P^. (123) 

At the critical field the (transverse) magnetization behaves as 

m z -m z c ~ (Q- fi c ) (ln|fi- fi c | - 1) (124) 

and the static (transverse) susceptibility \ z exhibits a logarithmic singu- 
larity (compare with Eqs. I|98l) . I|99|) L The bosonization approach results^ 
show that the same picture is valid for an arbitrary I z ^ provided 5 and 
7 are suitably renormalized. The results of bosonization approach elabo- 
rated for small 7 and S may be compared with numerical findings within 
non-perturbative regimes!^ From this analysis (0 < I z < 1) it was found 
that the critical line for 51 = is given by 

(5-7 (125) 

in accord with 1|123[) . For the behavior of the magnetization curves near the 
critical fields il c (I z = 1) a fair regime l|124|) was obtained. These outcomes 
suggest that some basic features of the fully interacting system (|122|) are 
captured by the free fermion picture. 
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